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Abstract 

miRNA levels are altered in pancreatic ductal adenocarcinoma (PDA), the most common and lethal pancreatic malignancy, 
and intact miRNA processing is essential for lineage specification during pancreatic development. However, the role of 
miRNA processing in PDA has not been explored. Here we study the role of miRNA biogenesis in PDA development by 
deleting the miRNA processing enzyme Dicer in a PDA mouse model driven by oncogenic Kras. We find that loss of Dicer 
accelerates Kras driven acinar dedifferentiation and acinar to ductal metaplasia (ADM), a process that has been shown to 
precede and promote the specification of PDA precursors. However, unconstrained ADM also displays high levels of 
apoptosis. Dicer loss does not accelerate development of Kras driven PDA precursors or PDA, but surprisingly, we observe 
that mouse PDA can develop without Dicer, although at the expense of proliferative capacity. Our data suggest that intact 
miRNA processing is involved in both constraining pro-tumorigenic changes in pancreatic differentiation as well as 
maintaining viability during PDA initiation. 
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Introduction 

Pancreatic ductal adenocarcinoma (PDA) appears to develop 
through a series of ductal precursor lesions, including the most 
common type, pancreatic intraepithelial neoplasia (PanIN). Both 
PanlNs and PDA exhibit Kras mutations, which may be an 
initiating event. Mouse models support this notion as targeted and 
persistent expression of mutant Kras in the mouse pancreatic 
epithelium both recapitulates the PanIN to PDA sequence 
observed in humans, and is required for disease maintenance 
[1,2,3]. Evidence from mouse models suggests that mutant Kras 
can contribute to PDA initiation by reprogramming acinar cells 
into a duct like lineage capable of becoming PanlNs via a process 
termed acinar to ductal metaplasia (ADM) [4,5,6]. This ability to 
change pancreatic plasticity may be an important step in PDA 
initiation, as acinar cells have been shown to be dramatically more 
sensitive to Kras dependent PanIN development compared to duct 
cells [7]. Since attempts at direct inhibition of oncogenic Kras 
have been generally unsuccessful [8], defining critical mediators of 
pro-tumorigenic differentiation and viability downstream of 
mutant Kras may represent alternative therapeutic approaches. 

MicroRNAs (miRNAs) are a class of small, non-coding RNAs 
that regulate gene expression post-transcriptionally. Aberrant 
miRNA levels are associated with tumor development, and lead 
to inappropriate expression of oncogenes and tumor suppres- 
sors [9]. Cancer associated miRNA levels result from the mis- 



expression of specific miRNAs as well as deregulation of the 
miRNA biogenesis pathway [10]. miRNAs are generally downreg- 
ulated in human cancers [11], and inhibited miRNA processing 
promotes tumorigenesis[12]. However, loss of miRNA processing 
may be incompatible with the development of some tu- 
mors [13, 14, 15], suggesting that critical thresholds of miRNA 
processing may be involved in maintaining viability during 
transformation. While miRNAs are misexpressed in PanlNs and 
PDA, and are known to regulate pathways that contribute to PDA 
initiation and progression[16,17,18,19], the role that miRNA 
processing plays in PDA development has not been investigated. 

By deleting the miRNA-processing enzyme Dicer in a Kras 
driven mouse model of PDA, we find that miRNA processing 
regulates both differentiation and viability during Kras driven 
pancreatic transformation. Dicer deletion promotes Kras driven 
loss of acinar identity and ADM but also results in increased levels 
of apoptosis and decreased expression of genes implicated in 
maintaining viability during PanIN development and in PDA. 
Surprisingly, we find that mouse PDA can develop in the absence 
of Dicer, although with decreased proliferation. Taken together, 
this work suggests a critical role for miRNA processing in Kras 
driven PDA initiation. 
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Results 

Dicer Loss Accelerates Kras Driven Ductal Metaplasia 

Dicer deletion in early pancreatic progenitors results in 
pancreatic agenesis and post-natal death [20]. Therefore, we 
tested if a distinct Pdxl driven Cre strain, Pdxl -Cre 1 "'' permitted 
pancreatic development in the setting of loss of Dicer function. 
Pdxl-Cre 1 "'" results in delayed developmental activity compared to 
the Pdxl -Cre driver employed by Lynn et al, resulting in 
recombination in a more restricted set of adult cells, specifically 
most acini, some endocrine cells, and rarely in duct cells [21]. 
Suggesting temporal requirements for Dicer in pancreatic 
development, Pdxl-Cre^"; DiceJ 1 "'^" 1 {Dice"'""") mice thrived and 
displayed grossly normal pancreatic development at pO even in the 
context of significandy decreased Dicer expression compared to 
control Pdxl-Cre 1 "''; Dicer"° x/+ mice (Dicer"' 1 ) (Figure SI A,B,E). 

At 3 weeks of age, acini, ducts, and islets were recognizable in 
Dicer 11 """' mice (Figure 1A), although there were some exocrine 
areas with decreased eosin staining. Immunofluorescence revealed 
normal distribution of the acinar marker amylase, duct marker 
CK19, and (3-cell marker insulin (Figure IB). Exocrine morphol- 
ogy was disturbed though, as we often found disorganized acini 
with small, fragmented cells not detected in controls (Figure IB, 
Figure S2), similar to disorganized acinar cells observed in somatic 
Dicer hypomorphs[22]. Pancreatic mass in Dicer H " m " mice was also 
reduced (Figure ID) in the setting of decreased Dicer expression 
(Figure 1C). Despite the changes in morphology, we did not 
observe cells co-expressing amylase and elevated levels of CK19 as 
observed in acinar cells undergoing Kras driven ADM (Figure S3), 
suggesting that the disorganized acinar cells were not actively 
undergoing ADM. 

To test the effect of Dicer loss on Kras driven pancreatic 
transformation, we generated Pdxl -Cre 1 *"; LSL-Kras G12D ; Dicer""' 
J 1 ™ (Kras; Dicer 11 '"") mice. Like other models targeting mutant Kras 
to the embryonic pancreas, Pdxl-Cre 1 " 1 '; LSL-Kras GI2D mice 
gradually develop ADM as well as PanlNs, with PDA arising in 
some animals after long latency[23]. Similar to such modelsfl], at 
3 weeks of age, pancreas mass in Pdxl-Cre 1 "'"; LSL-Kras G12D ; 
Dice/" x/+ (Kras; Dicer"") mice was slightly, but significantly, 
increased compared to Dicer 11 '' mice (Figure ID). The exocrine 
compartment in Kras; Dicer 11 '' mice appeared grossly normal, 
predominantly composed of amylase positive acinar cells, with rare 
areas of amylase negative ADM and low grade PanlNs expressing 
moderate to high levels of CK19 and Sox9 (Figure IF), ductal/ 
pancreatic embryonic progenitor markers characteristic of Kras 
driven ADM and PanIN formation [6] . In contrast, Kras; Dicer 11 """ 
pancreata displayed decreased Dicer expression (Figure 1C), were 
atrophic (Figure ID), and possessed widespread replacement of the 
exocrine compartment with metaplastic duct structures resembling 
rare ADM in Kras; Dicer 11 '' mice (Figure IE). Ductal metaplasia in 
Kras; Dicer 11 """' mice displayed decreased expression of amylase, low 
to moderate expression of CK19, and strong expression of Sox9 
(Figure IF). Sox9 accumulation was also observed in structures 
retaining some acinar morphology (Figure IF). Ductal metaplasia 
was not observed at pO in Kras; Dicer 11 "™" mice (Figure SID), 
suggesting that Dicer deficiency in the context of mutant Kras 
does not block embryonic acinar development, and ADM occurs 
between birth and 3 weeks of age. Therefore, Dicer loss in the 
context of mutant Kras dramatically accelerates ADM, a process 
that has been shown to both precede and promote Kras driven 
PanIN formation [4,6,24]. 



Dicer Loss Compromises Acinar Identity and Promotes 
Kras Driven Acinar to Ductal Reprogramming 

To explore why Kras; Dicer 11 "" 1 " mice undergo accelerated Kras 
driven ADM, we asked if Dicer deficient acinar cells inappropri- 
ately display properties associated with ADM, such as markers of 
acinar stress and loss of acinar identity. To focus on cells that had 
undergone Cre recombination we included a Cre inducible YFP 
allele expressed conditionally from the Rosa26 locus (R26-EYFP). 
YFP staining revealed that both the exocrine compartment in 
Pdxl-Cre 1 *"; DiceA'^'; R26-EYFP (Dicer""""'; YFP] mice and ADM 
in Pdxl-Cre 1 *"; LSL-Kras G12D ; Dicef" x,Jl ° x ; R26-EYFP {Kras; Dicer- 
Homo. Yppj m i ce derived from cells in which Cre was active, and not 
from expansion of an un-recombined population (Figure 2B, D). 

To determine if Dicer loss led to acinar stress associated with 
ADM, we examined expression of clusterin, a marker of stressed, 
de-differentiated acini[6,25]. Clusterin was mainly absent in YFP+ 
acinar cells of control Dicer 11 "; YFP and Kras; Dicer 11 "; YFP mice, but 
was present in rare YFP+ ADM in Kras; Dicer 11 "; YFP animals 
(Figure 2A, C). Suggesting that Dicer deficient cells are under 
stress observed in ADM, clusterin was not only uniformly observed 
in YFP+ ADM in Kras; Dicer 11 '"""; YFP mice, but also widely 
expressed in YFP+ acinar cells in Dicer 11 "" 1 "; YFP animals (Figure 2B, 
D). 

To test if acinar stress corresponded with decreased acinar 
identity, we modified a sorting protocol previously developed for 
isolation of pancreatic progenitors [26] to specifically separate 
differentiated acinar and ductal cells. We found that in the adult 
pancreas, adult acini and ducts express the epithelial marker 
CD49f, while differentiated ducts also express CD 133 (Figure S4B, 
C). RT-PCR for Amylase and Mistl, and CK19 and Sox9, enriched 
in acinar and duct cells, respectively, revealed that sorting based 
on CD49f and CD 133 allowed for efficient separation of these 
populations (Figure S4A, 4D). Expression of Notch effectors Hesl, 
Heyl, and Hey2, restricted to centroacinar and terminal duct cells 
in the adult pancreas [2 7], appeared to segregate with the CD49f+ 
CD133+ duct population (Figure S4E). 

YFP+, CD49f+CD133- cells were sorted from all 4 genotypes. 
At 3 weeks of age Dicer levels and expression of 2 abundandy 
expressed miRNAs were considerably reduced in cells sorted from 
Dicer 11 '"""; YFP and Kras; Dicer 11 "" 1 "; YFP mice compared to controls 
with and without Kras (Figure 2E left). To determine the effect of 
Dicer loss on exocrine differentiation, we analyzed expression of 
the acinar markers Amylase and Mistl and the duct markers 
CK19 and Sox9, in YFP+, CD49f+CD133- cells. Sorted cells from 
Dicer 11 """; YFP mice displayed considerably decreased expression of 
Amylase and modestly decreased Mistl compared to cells from 
Dicer "; YFP mice (Figure 2E middle). Although expression of 
markers of acinar differentiation were reduced, we did not observe 
a consistent increase in the ductal markers CK19 or Sox9 in these 
cells (Figure 2E right), similar to the lack of widespread CK19 
misexpression noted in Figure ID. In CD49f+CD133- cells from 
Kras; Dicer"'"""; YFP mice however, Amylase and Mistl expression 
was further decreased, and CK19 and Sox9 expression was 
increased, compared to cells from all other genotypes, consistent 
with histological increase in ADM. 

Interestingly, we also observed a modest, but not statistically 
significant, trend toward reduction in Mistl expression in cells 
sorted from Kras; Dicer""; YFP mice compared to Dicer""; YFP 
controls, suggesting that mutant Kras alone may begin compro- 
mising acinar identity before the induction of ductal genes or 
morphology. Taken together, this data suggests that Dicer loss 
compromises acinar identity and induces a stressed state that 
accelerates Kras driven ductal reprogramming and ADM. 
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Figure 1 . Dicer loss accelerates Kras driven ADM. (A). Hematoxylin and Eosin (H&E) staining of 3 week old pancreas from Dicer"" and Dicer" 0 " 10 
mice. (A) acini, (D) duct, and (I) islets. Scale bar 100 ^M. (B) Immunofluorescence for amylase (red), CK19 (green), and insulin (blue) in 3 week old 
Dicer" et and Dicer" 0 '" 0 mice. Scale bar 100 ^M. (C). Reduced Dicer expression at 3 weeks in RNA extracted from pancreas tissue of indicated 
genotypes. Mean ± SD, n = 3. (D). Pancreas mass: body weight ratio at 3 weeks in the indicated genotypes. P-values are calculated from two-tailed, 
unpaired t-tests comparing Dicer""' and indicated genotypes. (Mean ± SD, Dicer" e ' n = 1 7, Dicer" omo n = 14, Kras; Dicer"" n = 9, Kras; Dicer" 0 "' 0 n = 1 4). 
(E). H&E staining of 3 week old Kras; Dicer""' and Kras; Dicer" 0 " 10 mice. Insets: E-left, Rare focus of metaplasia and PanIN; right, ductal metaplasia in 
Kras; Dicer" 0 " 10 mice. Scale bar 1 00 \iM. (F) Immunohistochemistry for Amylase, Sox9, CK1 9 in 3 weeks old Kras; Dicer""' and Kras; Dicer" omo mice. Scale 
bar 100 u.M. 

doi:1 0.1 371 /journal.pone.0095486.g001 



Dicer Loss does not Accelerate PanIN or PDA 
Development 

Loss of acinar identity and ADM accelerates PanIN develop- 
ment [6,24,28,29] and accelerated ADM and PanIN development, 
induced by acute and chronic pancreatitis for example, can result 
in decreased PDA latency [30,31,32]. Therefore, we expected to 
observe faster PanIN and PDA development in Kras; Dicer 1 """ mice 



compared to control Kras; Dicer '' mice. Despite the dramatic 
difference in ADM observed at 3 weeks in Kras; Dice/ 1 """' 
compared to Kras; Dicer 11 '' mice, quantification of PanIN lesions 
in a small cohort of mice at 9 weeks revealed no statistically 
significant difference in the frequency of Alcian Blue positive 
ductal lesions, suggesting that Dicer loss does not accelerate PanIN 
development (Figure 3A, B, C). Long-term disease progression also 
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Figure 2. Dicer loss promotes Kras driven loss of acinar identity and ductal reprogramming. (A-D) Clusterin (blue) and YFP (green) 
staining in 3 week old Dicer""; YFP (A), Dicer""™; YFP (B), Kras; Dicer"*'; YFP (C), and Kras; Dicer" 0 ™; YFP mice (D). Scale bar 1 00 \M (E) RT-PCR analysis 
of Dicer, miRNA-141 and 216b (left), acinar enriched genes Amylase and Mistl (center), and duct enriched genes CK19 and Sox9 (right) of YFP+, 
CD49f+, CD133- cells from Dicer""; YFP (n = 4), Dicer""""'; YFP (n=4), Kras; Dicer" a ; YFP (n = 3), and Kras; Dicer""™; YFP (n = 3) mice. Mean ± SD. P- 
values are calculated from two-tailed, unpaired t-tests comparing Dicer""' and Dicer"° m CK19 and Sox9 expression values. 
doi:1 0.1 371 /journal.pone.0095486.g002 
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converged independent of conditional Dicer status. We found no 
statistically significant difference in survival of a small cohort of 
aged Kras; Dker Hmm and Kras; Dicer 11 '' mice (Median survival 336 
versus 441.5 days, / 2 = .2271, p = 0.6637, calculated using the 
Logrank test) (Figure 3F). 6/7 of Kras; Dice/ 1 ' 1 mice developed 
locally or widely invasive PDA that was generally moderately 
differentiated, but included one undifferentiated cancer. 4 out of 6 
Kras; Dice/ 1 "" 1 " mice developed moderately differentiated, locally or 
widely invasive, PDA with similar frequency of metastasis and 
cellular morphology (Figure 3D, E, Table SI). 

Dicer deleted cells have been shown to be outcompeted by 
unrecombined cells in other endoderm derived organs, in the liver 
for example [33], and elimination of recombined cells could 
contribute to the lack of expected difference in disease progression 
in Kras; Dicer 110 ""' versus Kras; Dicer 11 '' mice. In contrast to the 3 
week time point when normal acinar cells were much less frequent 
m Kras; Dicer 110 ™ versus Kras; Dicer 11 " mice, considerable areas of 
normal acinar tissue were found in Kras; Dicer 11 " 1 ™ mice at this 
time point (Figure 3 A, B). We also noted the presence of fat 
replacing pancreas tissue in some Kras; Dice/ 1 "'" 0 animals 
(Figure 3B). Examination of tissue YFP distribution in the 
pancreas of Kras; Dicer 11 " 1 ; YFP versus Kras; Dicer 110 " 1 "; YFP animals 
revealed a notable decrease in YFP distribution at 9 weeks (Figure 
S5), suggesting repopulation with cells that had not undergone 
recombination. Therefore, despite initially promoting ADM that 
normally accelerates PanIN and PDA development, PDA 
progression is not substantially enhanced in the absence of Dicer 
and may be associated with selection against Dicer deficiency. 

Dicer Loss Compromises Viability during Kras Driven 
ADM 

Dicer loss leads to increased apoptosis in a number of 
endodermal organs (e.g. liver and intestine [33,34]. To address if 
the lack of predicted acceleration of PanIN development in Kras; 
Dicer 11 ""™ mice might involve increased cell death, we audited 
TUNEL+/YFP+ cells at 3 weeks of age. In contrast to the rare 
double positive cells in Dicer 11 ''; YFP and Kras; Dicer 11 " 1 ; YFP mice, 
Dicer; YFP mice displayed a higher relative rate of TUNEL+/ 



YFP+ cells (Figure 4B, E). Furthermore, we found an additional 
increase in double positive cells in Kras; Dicer 11 "" 1 "; YFP mice 
compared to Dicer""; YFP and Kras; Dicer 11 "; YFP mice(Figure 4D, 
E). Interestingly, Dicer loss appeared to synergize with mutant 
Kras to promote cell death, as Kras; Dice/ 10 " 10 ; YFP mice exhibited 
significandy more apoptosis than Dicer 110 " 10 ; YFP cells (Figure 4E). 
Therefore, while intact Dicer processing constrains Kras driven 
acinar to ductal reprogramming, some Dicer dependent signals 
may be required for maintaining viability during Kras driven 
metaplasia. 

We performed expression analysis to screen for Dicer dependent 
factors potentially involved in maintaining viability during Kras 
driven ADM. To better synchronize ADM in Kras; Dicer 11 " and 
Kras; Dicer 110 " 10 mice, we utilized caerulein induced pancreatitis. As 
predicted from previous studies [6] , caerulein treatment in control, 
3 weeks old Kras; Dicer 11 " mice led to the development of duct like 
structures 2 days after treatment (Figure S6A), and widespread 
replacement of the exocrine compartment with ADM, PanlNs, 
and fibrosis, 21 days after treatment (Figure S6C). Mirroring the 
lack of accelerated PanIN development between 3 and 9 weeks in 
Kras; Dicer""'"" mice (Figure 3F), we observed less metaplasia and 
considerably more normal acinar tissue 21 days following 
caerulein treatment in Kras; Dicer 11 "" 10 mice (Figure S6D). RNA- 
Seq Analysis of pooled RNA from FACS isolated YFP+, CD49f+, 
CD 133- acinar cells from 3 weeks old Kras; Dice/ 1 "; YFP (3 mice) 
and Kras; Dicer ; YFP (5 mice) mice 2 days after caerulein (Figure 
S6E) revealed — 1 20 genes significantly downregulated between 
mutant and control. Supporting the observation that loss of Dicer 
compromises acinar differentiation, 42 genes associated with 
terminally differentiated acinar cells were significantly downreg- 
ulated (Figure 4F, Table S2, complete data set Table S3). Also 
reduced were genes implicated in maintaining viability during 
regeneration or tumorigenesis, including Nuprl/p8, Agr2, Itih4, 
and members of the Reg3 protein family (Figure 4F, Table S2). 
Recently, both Nuprl and Agr2 have been shown to be 
overexpressed during PanlN-PDA progression and to play a role 
in PanIN development and maintaining PDA viability [35, 36,37]. 
Reg3 family protein expression is activated in response to 
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Figure 3. Kras driven PanIN and PDA progression in Dicer conditional mice. (A-B) Representative Alcian Blue staining at 9 weeks in Kras; 
Dicer"" (A) and Kras; Dicer" omo mice (B). f: fat. Scale bar 1 00 uM. C, Quantification of alcian blue positive lesions of 9 weeks old Kras; Dicer""' and Kras; 
Dicer" om ° mice. Lines indicate means. P-value calculated from a two-tailed, unpaired t-test. (D,E) Representative histology of PDA from Kras; 
Dicer" e '(D) and Kras; Dicer" 01 " 0 mice (E). Scale bar 1 00 uM. (F) Survival curve of Kras; Dicer"" (n = 7) and Kras; Dicer Homo (n = 6) mice. 
doi:1 0.1 371 /journal. pone.0095486.g003 
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Figure 4. Dicer loss compromises viability during unconstrained Kras driven ADM. (A-D) TUNEL (blue) and YFP (green) staining in 3 
weeks old Dicer"*'; YFP (A), Dicer" 0 ™; YFP (B), Kras; Dicer"*; YFP (C), and Kras; Dicer" 0 ""'; YFP mice (D). Scale bar 100 nM. E. Quantification of TUNEL+ 
YFP+ cells in (A-D). Mean ± SD. Dicer""'; YFP (n = 3), Dicer" omo ; YFP (n = 4), Kras; Dicer""'; YFP (n = 3), and Kras; Dicer" omo ; YFP mice (n = 4). P-value 
calculated from a two-tailed, unpaired t-test of TUNEL+YFP+ cells in Dicer" omo ; YFP and Kras; Dicer" omo ; YFP mice. F. Acinar enriched and viability 
associated genes from expression analysis (left). RT-PCR analysis of Nuprl, Agr2, Itih4, and Reg3b in acinar enriched cells from Dicer"*'; YFP (n=4), 
Dicer" omo ; YFP (n = 4), Kras; Dicer"*'; YFP (n = 3), and Kras; Dicer" orno ; YFP (n = 3) mice (right). Mean ± SD. P-value calculated from a two-tailed, unpaired 
t-test comparing Agr2 expression in sorted cells from Dicer"*'; YFP and Kras; Dicer"* 1 ; YFP mice. 
doi:10.1371/journal.pone.0095486.g004 



pancreatitis [38], and family member Reg3b (also known as Papl) 
specifically is implicated in maintaining acinar viability and can 
act downstream of Nuprl [39,40]. Although Itih4 has not been 
implicated in PDA development, it has been shown to be involved 
in maintaining viability downstream of the critical PDA mediator 
c-myc[41] in a model of liver cancer[42]. RT-PCR analysis of 
these candidate genes in sorted acini from 3 week-old mice 
revealed that their expression inversely correlated with apoptosis 
(Figure 4F, right). Nuprl, Agr2, Itih4, and Reg3b expression was 
reduced in sorted cells from Kras; Dicer H ° m °; YFP mice compared to 
control Dicer Hlrt ; YFP and Kras; Dice/ 1 "'; YFP animals. Also, Agr2, 
Itih4, and Reg3b expression was reduced in sorted cells from 
Dicer"'""''; YFP mice compared to controls. Interestingly, Agr2 was 



also significantly increased in cells sorted from Kras; Dicer Het ; YFP 
mice compared to cells from Dicer 11 "'; YFP mice, suggesting that 
Kras may contribute to Agr2 upregulation even before ADM 
occurs. Therefore, Dicer expression is required for the mainte- 
nance of the mature acinar differentiation state and for expression 
of pro-viability genes during Kras driven ADM. 

Mouse PDA can Develop in the Absence of Dicer 

Dicer levels have been shown to be important regulators of 
tumorigenesis, acting as a haploinsufficient tumor suppressor in 
the context of mutant Kras in both lung tumors and sarcomas, as 
well as in retinoblastoma [1 3, 14] . However, tumors in these models 
appear to result from cells that have retained expression of at least 
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one unrecombined conditional allele. To determine if Dicer 
deficient cells were capable of contributing to Kras driven PDA we 
tested Dicer expression and genomic recombination from 3 PDA 
cell lines generated from Kras; Dicer 11 ' 1 and Kras; Dicer 11 "" 1 " mice. By 
RT-PCR we observed similar Dicer expression in all three PDA 
cell lines derived from Kras; Dicer 11 "' mice (Figure 5A). Due to the 
high rate of apoptosis observed at 3 weeks of age, the reduced 
expression of genes that support PanIN development, and the lack 
of accelerated disease progression, we expected to observe 
retention of at least one allele in PDA lines derived from Kras; 
Dicer 11 """' PDA. Surprisingly, we found that one line displayed 
dramatically decreased Dicer expression, while the two remaining 
lines expressed either a similar, or higher, level as Kras; Dicer 11 ' 1 
derived lines. Allele specific PCR revealed that Dicer recombina- 
tion correlated directiy with expression (Figure 5B: "Undel", no 
conditional Dicer allele recombination; "Hemi", recombination of 
one allele; "Homo", recombination of both alleles). RT-PCR for 3 
mature miRNAs reported to be over-expressed in human 
PDA[43,44] revealed reduction in all 3 miRNAs in the 
homozygous deleted line to less than 95% of that expressed in 
the undeleted line (Figure 5C). Therefore, Kras driven PDA in the 
mouse appears able to evade negative selection during disease 
initiation and develop in the absence of Dicer. As predicted from 
models demonstrating that downregulation, but not complete loss 
of miRNA processing, can enhance tumorigenesis[12,13], the 
hemizygous deleted line grew fastest, while the homozygous 
deleted line slowest (Figure 5D). All 3 cell lines were also able to 
form tumors when implanted subcutaneously in immune deficient 
mice, with kinetics mirroring the cell culture growth rates 
(Figure 5E). Allele specific PCR revealed that tumors did not 
result from cells that underwent spontaneous recombination 
(UnDel, Hemi), or from a contaminating subpopulation of Dicer 
competent cells (Homo), although we did observe the presence of 
the WT Dicer amplicon in tumors derived from the Hemi and 
Homo lines, indicative of host stromal recruitment (Figure 5F). 
Therefore, although Dicer elimination impairs viability during the 
earliest stages of Kras driven neoplasia, loss of Dicer is not 
mutually exclusive with pancreatic transformation. 

Discussion 

Dicer function is essential for the development of several organs, 
including endoderm-derived tissues such as the liver and the 
gut[33,34]. Maintenance of Dicer expression in early Pdxl 
positive pancreatic progenitors is required for development of 
both the pancreatic exocrine and endocrine compartments [20]. 
Here, we find that Dicer function remains important throughout 
pancreatic development and plays a role in regulating acinar 
identity and viability. Our model utilizes a driver line that initiates 
Dicer recombination in the developing pancreas at a later 
embryonic stage than the Cre driver used by Lynn and 
colleagues [21]. We observe minimal gross effects on pancreatic 
development, suggesting stage sensitive requirements for miRNA 
processing on establishing and expanding exocrine and endocrine 
progenitors. However, the acinar cells that do develop are unstable 
both in regards to their terminal differentiation and viability. 
Elements of acinar identity are downregulated, and cell stress and 
apoptosis are both increased. Therefore, competent miRNA 
processing appears to remain important even in the mature 
exocrine compartment. 

Because Dicer deletion in this model permits pancreatic 
development we were able to explore the role of Dicer function 
in Kras mediated PDA development. Mouse models have revealed 
that oncogenic Kras can act as a "master regulator" of PDA 



development, establishing lineages that can give rise to PanlNs 
and PDA and remaining critical for progression [1,2, 3]. Consid- 
erable evidence suggests that acinar cells can give rise to PanlNs 
by undergoing ADM, a process during which acinar cells lose 
terminal differentiation at the expense of a de-differentiated, duct 
like state[4,5,6,7]. Kras driven ADM occurs gradually, but can be 
accelerated by compromising acinar differentiation. Insults such as 
pancreatitis, which causes regeneration associated dedifferentia- 
tion[6,45], activation of progenitor associated pathways (e.g. 
Notch[24]), and inactivating genes that maintain the acinar 
state [28,29,46] all dramatically accelerate Kras driven ADM and 
PanIN development. Therefore, loss of acinar identity may be a 
key barrier to Kras dependent specification of acinar derived PDA 
precursors. Our data suggest that loss of Dicer function removes 
the differentiation barrier for Kras driven ADM. ADM occurs at 3 
weeks in Kras; Dicer 11 ""'" mice, a time point when the consequences 
on exocrine differentiation and morphology due to mutant Kras 
alone are minimal. While Dicer loss compromises acinar 
differentiation, it does not appear singly sufficient to induce ductal 
reprogramming, indicating synergy between Kras and Dicer loss 
to drive ADM. Supporting this role for the miRNA pathway in 
maintaining acinar differentiation, Prevot and colleagues [47] 
have recently also shown that Dicer deficient acini display 
compromised expression of acinar genes and upregulated expres- 
sion of genes involved in liver development, including the 
transcription factor HNF6/onecutl. Indeed, we find HNF6/ 
Onecutl to be highly expressed in sorted acinar cells from 
caerulein treated Kras; Dicer 11 "" 1 " versus Kras; Dicer 11 '' mice (Table 
S3), suggesting that inappropriate HNF6 /Onecutl expression 
may also be involved in Kras driven ADM. 

Despite inducing unconstrained Kras driven ADM, Dicer loss 
neither accelerates PanIN nor PDA development. This disconnect 
is potentially due to increased apoptosis, and thus selection against 
Dicer deleted cells, demonstrating that maintenance of appropri- 
ate levels of miRNA biogenesis is required for viability during Kras 
driven ADM. These opposing roles, pro-tumorigenic de-differen- 
tiation on one hand, cell death on the other, suggest checkpoints 
might be in place to eliminate miRNA deficient cells undergoing 
potentially deleterious, transformation sensitive de-differentiation. 
Supporting this notion, a similar pattern of compromised 
differentiation in the setting of increased apoptosis has been 
observed in other Dicer loss of function models. In the liver, for 
example, Dicer loss leads to the inappropriate activation of pro- 
proliferative genes expressed in fetal liver development, accompa- 
nied by overwhelming levels of cell death [3 3]. Synthetic lethality 
between pro-oncogenic signaling, namely tumor suppressor loss, 
and Dicer function has also been observed in retinoblastoma (Rb), 
p53, Dicer null retinal progenitor cells [48]. This lethal interaction 
resembles the synergistic cell death we observe between loss of 
Dicer function and mutant Kras during unconstrained ADM. 
Increased cell death also corresponds to decreased expression of 
genes that have been experimentally validated to support viability 
during Kras driven PanIN development and in PDA cells, such as 
Nuprl and AGR2 [35,36,37,49]. If Dicer dependent regulation of 
these genes depends on specific miRNAs throughout PDA 
development, targeting those miRNAs directly could also be a 
route to undermining PDA viability. 

Dicer dependent cell death during Kras driven ADM suggests 
that pancreatic cells exhibit similar negative selection against Dicer 
loss similar to that seen in a number of tumor models, 
[13,14,15,48,50]. In line with this observation, 2 of 3 cell lines 
derived from PDA in Kras; Dicer 11 "" 1 " mice maintained at least one 
intact Dicer allele. Surprisingly though, one of the 3 cell lines 
derived from PDA in Kras; Dicer"'""" mice displayed recombination 
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Figure 5. Dicer loss can be tolerated during mouse PDA development. A. Dicer RT-PCR in PDA cell lines generated from Kras; Dicer"*' (blue) 
and Kras; Dicer" 0 "' 0 (red) mice. Each bar represents data from 3 independent wells of each cell line. Mean ± SD. B. PCR detecting WT, unrecombined 
(UNREC), and deleted (DEL) Dicer genomic locus in cell lines from A. Controls are amplified DNA from mouse embryonic fibroblasts (MEFs) of 
indicated Dicer genotype following adenoviral Cre treatment. C. miRNA QPCR for miRNA-16, 107, 191 in cell lines derived from Kras; Dicer" 0 " 10 mice. 
Each bar represents 3 independent wells of each cell line. Mean ± SD. D. Growth of PDA cell lines derived from Kras; Dicer""™" mice. Each point 
represents values from 3 independent wells. Mean ± SD. P-value calculated from a two-tailed, unpaired t-test comparing the number of cells at 72 
hours in cultures of Undel and Homo cells. E. Subcutaneous tumor growth of PDA cell lines derived from Kras; Dicer"""' 0 mice. Mean ± SD at each 
timepoint. Undel (n = 7), Hemi (n = 7), Homo (n = 1 2). F. PCR detecting WT, UNREC, and DEL Dicer genomic locus in subcutaneous tumors at time of 
sacrifice in E. Controls are Dicer flox/+ MEFs treated with adenoviral Cre as in B. Tumor DNAs are preceded by amplification of DNA from parental cell 
lines. 

doi:1 0.1 371 /journal.pone.0095486.g005 



of both Dicer alleles and appears to have resulted from a tumor 
that developed in the absence of Dicer. Processing of some 
miRNAs has been shown to occur in a Dicer independent fashion, 
including miR-451 [51,52], a miRNA that has been noted to be 
abundant in PDA fine needle aspirate samples [5 3]. Although 
beyond the scope of this study, compensation by Dicer indepen- 
dent miRNAs could possibly be involved in maintaining viability 
in the absence of Dicer function. The invasive histology of PDA 
from Kras; Dicer ° m ° mice appeared similar, transformation without 
Dicer occurs at the expense of some pro-tumorigenic cell biology, 
as observed in clonally selected Dicer null sarcomas [50]. Cells 
possessing one Dicer allele and intermediate Dicer expression 
proliferated faster, and grew larger tumors when transplanted into 
immune deficient mice, compared to both Dicer deficient cells and 
cells retaining both Dicer copies. Therefore, PDA, and other 
mutant Kras driven tumors [12], may depend on maintaining 



miRNA processing above a critical threshold required for viability 
and proliferation. While monoallelic Dicer deletions have been 
reported[13], and variable Dicer levels have recently been 
observed in PDA cell lines [54], extensive characterization of 
Dicer and components of the miRNA biogenesis machinery in 
PDA has not been performed. 

In conclusion, we establish intact miRNA processing as an 
important mediator of Kras driven PDA development. We find 
that intact Dicer function constrains Kras driven ADM by 
maintaining acinar differentiation but is also required for cell 
viability during Kras driven metaplasia. Surprisingly though, 
mouse PDA can develop without Dicer. Future work understand- 
ing the relationship between thresholds of miRNA processing and 
differentiation and viability could help better define targets for 
PDA therapy. 
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Materials and Methods 

Mouse Lines 

Experimental animals were generated by intercrossing mice 
bearing Pdx-Cre^"' [21], LSL-Kras G12D [1], Dice/'" [33] and R26R- 
EYFP[6] alleles. All mouse experiments were performed under the 
approval of the University of California, San Francisco Institu- 
tional Care and Use of Animals Committee (IACUC, Institutional 
PHS Assurance #A3400-01), and all efforts were made to 
minimize animal suffering. 

Immunohistochemistry and Immunofluorescence 

Pancreata were fixed overnight in zinc-containing neutral- 
buffered formalin (Anatech LTD), embedded in paraffin, cut into 
5-um-thick sections, and mounted on Superfrost Plus slides (Fisher 
Scientific). Sections were subjected to hematoxylin and eosin 
(H&E), immunohistochemical, and immunofluorescent staining as 
described [6] . The following primary antibodies were used: rabbit 
anti-amylase (1:300; Sigma), rat anti-CK19 (TROMAIII, 1:200 
dilution; developed by Dr. Rolf Kemler [Max-Planck Institute of 
Immunobiology, Freiburg, Germany] and obtained from the 
Hybridoma Bank at the University of Iowa), goat anti-clusterin 
(1:200; Santa Cruz), rabbit anti-Sox9 (1:1000; Chemicon), chicken 
anti-GFP (1:200; Abeam), biotinylated anti-Cd49f (1:50 
eBioscience), and biotinylated anti-CD 133 (1:50, eBioscience). 
For immunohistochemistry, biotinylated secondary antibodies 
were used at a 1:200 dilution and peroxidase conjugation was 
performed with the streptavidin based ABC kit (V ector labs). 3-3- 
Diaminobenzidine tetrahydrochloride (V ector Labs) was used as a 
chromogen. Bright-field images were acquired using a Zeiss Axio 
Imager Dl scope. For immunofluorescence detection of primary 
antibodies, appropriate Alexafluor conjugated secondary antibod- 
ies were used at a 1:200 dilution. Confocal images were collected 
on a Leica SP5 microscope. Alcian blue staining was performed as 
described [6], and quantified by collecting 100X images encom- 
passing at least 1 complete tissue section from each mouse and 
scoring the number of Alcian Blue positive lesions per mm 2 of total 
pancreatic area (measured by outlining pancreatic tissue with the 
Axiovision software package (Zeiss)). 

TUNEL Assay 

TUNEL staining was performed with the Apoptag Fluorescein 
in situ Apoptosis Detection Kit (Chemicon) according to manu- 
facturer's instructions. The number of TUNEL+ cells per unit 
YFP+ area were quantified by imaging at least 1 entire tissue 
section at 100X on an InCell image automated microscope (GE). 

Cell Sorting 

Single cell suspensions of 3 week old pancreata were generated 
with a modified version of a protocol established by Sugiyama and 
colleagues, 2007 [26]. Briefly, pancreatic lymph nodes were 
removed, and pancreata were minced and sequentially incubated 
with collagenase D (1 mg/ml in Hanks buffered saline solution, 
Roche), trypsin (0.05%, Invitrogen), and dispase (2 U/ml, 
Invitrogen). Dnasel (100 ug/ml, Sigma) was added during all 
enzyme incubations. Cells were washed with PBS between the 
collagenase and trypsin steps, and with FACs buffer (2% FBS, 
10 mM EGTA, in PBS) between the trypsin and dispase steps. 
Suspensions were then filtered through 40 uM mesh, subjected to 
FC block (1:200), and incubated with PE-conjugated anti-Cd49f 
(1:50, Ebioscience), PE-CY7-conjugated anti-CD45 (1:500, 
Ebioscience), and biotin conjugated anti-CD 133 antibodies 
(1:50, Ebioscience). Streptavidin conjugated to APC was then 
applied to detect CD 133 positive cells (1:100, Ebisocience). Cells 



were washed with FACs buffer between block, primary and 
secondary antibody incubations. DAPI (300 nM) was used as a live 
cell marker. Cell sorting was performed on a FacsARIAII (Becton 
Dickinson). 

RT-PCR 

RNA was extracted from sorted cells using the RNAqueous RNA 
extraction kit (Ambion) according to modified manufacturer's 
instructions to include small RNA fractions, and from tissue with 
the RNEasy kit (Qiagen) as described [6]. cDNA was synthesized 
using Superscript II Reverse Transcriptase (Invitrogen). Taqman 
RT-PCR was performed using inventoried probes for mouse Amylase 
2 (Mm02342486_mH), Mistl (Mm00487695_ml), CK19 
(Mm00492980_ml), Sox9 (Mm00448840_ml), Hesl (Mm0134- 
2805_ml), Heyl (Mm00468865_ml), Hey2 (Mm00469280_ml), 
HeyL (Mm00516555_mL) (Applied Biosystems), Nuprl (Mm0049- 
8104_ml), Reg3b (Mm00440616_gl), ITIH4 Mm00497648_ml), 
and Agr2 (Mm005()7853_ml). Expression levels were normalized 
using a custom primer/ probe set for Cyclophilin (generously provided 
by the Genome Analysis Core at the UCSF Helen Diller Family 
Comprehensive Cancer Center: Primer 1: GGCCGATGAC- 
GAGCCC, Primer 2: TGTCTTTGGAACTTTGTCTGCAA, 
Probe: TGGGCCGCGTCTCCTTCGA). RT-PCR for Dicer was 
performed using SYBR GREEN master mix (Applied Biosystems) 
with primers specifically designed to detect deletion of the conditional 
Dicer allele [33]. Dicer expression was normalized to Gus levels 
detected by RT-PCR with SYBR GREEN using primers designed 
by Heiser and colleagues (2008) [21]. miRNA RT-PCR was 
performed with the Taqman miRNA reverse transcription kit 
(Applied Bisosystems) and inventoried assays for mouse miR-141 
(000463), miR-216a 002220), miR-16 (000391), miR-107 (000443), 
and miR-191 (002299). miRNA expression was normalized to 
SnoRNA-202 levels (Applied Biosystems, 001232) 

Caerulein Treatment 

PdxCre^"; LSL-Kras G12D ; Dice/ 1 "''* and PdxCre^"; LSL-Kras G,JD ; 
Dicef°' /Jl °'\ mice (both with and without the R26-EYFP allele) 
were treated with caerulein as described [6] , with 2 sets of 6 hourly 
i.p. caerulein injections (American Peptide Company, 50 ug/kg) 
on alternating days separated by 24 hours. The final day of 
caerulein injection was considered day 0. 

RNA Deep Sequencing 

500 ng of total RNA isolated from sorted pancreatic cells was 
used to generate libraries with the Illumina TruSeq RNA sample 
preparation kit. Each library was diluted to approximately 10 pM 
prior to loading and sequenced using a Hi-seq 2000 instrument 
generating 50 base pair reads. Library preparation and sequencing 
was performed by the Duke IGSP Genome Sequencing and 
Analysis Core. The 50 bp long reads were mapped to the mouse 
genome (NCBI37/mm9) and exon-exon junctions using TopHat 
[55] version 1.3.1 with default parameters. Genome information 
was downloaded from the UCSC Genome Browser [56]. We 
filtered TopHat alignments and kept reads (and their locations) 
that could be mapped uniquely. Per gene, we extracted read 
counts by using samtools [57] version 0.1.13. Read counts were 
normalized to the length of the gene and the total number of 
mapped reads to obtain a standardized count (RPKM, reads per 
kilobase of exon model per million mapped reads). Due to the lack 
of replicates, a Fisher's Exact Test was used to search for 
differentially expressed genes. P-values were calculated by 
comparing the normalized expression level of each gene to the 
normalized expression level of two housekeeping genes (Gusb, 
Rpl23). Complete data set is included as Table S3. 



PLOS ONE | www.plosone.org 



9 



May 2014 | Volume 9 | Issue 5 | e95486 



miRNA Processing and Pancreatic Cancer Development 



Establishment of PDA Cell Lines and Subcutaneous 
Tumor Growth Assays 

PDA bearing mice we sacrificed and fibrotic tumor areas were 
removed, minced with scissors, and sequentially digested with 
collagenase D, trypsin, and dispase as described in the cell sorting 
methods. Suspensions were filtered through 40 uM mesh and both 
trapped chunks and flow through were plated on collagen I coated 
plates (BD). Cells were passed until morphologically clear from 
fibroblasts and maintained on collagen I coated plates in DMEM 
containing 10%FBS and penicillin and streptomycin. For subcu- 
taneous tumor growth assays, 5x10 tumor cells resuspended in 
PBS were mixed 1 : 1 with growth factor reduced matrigel (BD) and 
injected subcutaneously on the back of 6 week old immunodefi- 
cient NOD/SCID/gamma null mice. Tumor volume was 
estimated with calipers using the formula 0.5(Tength*(Width) 2 ). 
Tumors were removed when they reached ~500 mm 3 . 

Supporting Information 

Figure SI Deletion of Dicer with Pdxl-Cre Late permits 
pancreatic development. (A-D) Grossly normal pancreas 
histology in Dicer Het (A), Dice/ 1 ""™ (B), Kras; Dicer"" (C), and Kras; 
Dicer"'""" (D) mice at pO. Scale bar 100 uM. (C, D). RT-PCR 
reveals efficient pancreatic Dicer deletion in Dicer""™" versus 
Dicer"" mice at pO. Mean ± SD. n = 3. E. Reduced Dicer 
expression at pO in RNA extracted from Dicer""' and Dicer" 0 " 10 
pancreas. Mean ± SD, n = 3. 
(TIF) 

Figure S2 Acinar morphology in 3 weeks old Dicer 11 "* 
and Dicer 11 """' pancreas. (A) Normal acinar morphology and 
amylase distribution in 3 weeks old Dicer"'' mice. (B) Disturbed 
acinar morphology and fragmented amylase positive structures 
(arrowheads) in 3 weeks old Dicer""""' mice. Images are magnified 
regions from Figure IB. Note the absence of CK19 expression in 
Dicer"'""" acinar cells. 
(TIF) 

Figure S3 Loss of amylase and increased CK19 expres- 
sion in mouse ADM. (A-D) Immunofluorescence staining of 
YFP, Amylase, and CK19 in a 3-weeks old Kras; Dicer""; YFP 
mouse. Arrowhead indicates a YFP+ structure developing ductal 
morphology with a region of low amylase and high CK19 
expression. 
(TIF) 

Figure S4 FACS enrichment of pancreatic acinar and 
ductal compartments. (A) Representative flow cytometry plot 
of a pancreas from an adult, 6 week old mouse, dissociated into a 
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single cell suspension. Viable, non-hematopoietic cells (DAPI- 
CD45-, left panel) were gated and further analyzed for expression 
of CD49f and CD133 (right panel). (B, C) CD49f and CD133 
staining in 6 week old mice. Scale bar 100 uM. (D) RT-PCR 
analysis of acinar and ductal markers from CD49f+CD133- and 
double positive CD49f+CD133+ cells. (E) Analysis of Notch 
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(TIF) 

Figure S5 Loss of YFP positive cells at 9 weeks in Kras; 

Dicer"'"""; YFP mice. (A) Widespread YFP expression in 9 
weeks old Dicer""; Kras; YFP mice. (B) Reduced and heteroge- 
neous YFP expression in a 9 weeks old Dicer""""'; Kras; UP mouse. 
(TIF) 

Figure S6 Isolating acinar cells for determination of 
Dicer dependent genes during Kras driven ADM. (A,B). 

H&E staining of Kras; Dicer"" (A), and Kras; Dicer""""' (B) pancreas 2 
days after caerulein treatment. Scale bars 100 mM (C,D). H&E 
staining of Kras; Dicer"" (C), and Kras; Dicer"""" 1 (D) pancreas 2 1 
days after caerulein treatment. Scale bars 100 uM (E). 
Representative sorting profiles of Kras; Dicer""; EYFP and Kras; 
DicerHomo; YFP 2 days after caerulein. Gated CD49f+, CD 133- 
populations in right panels were collected for analysis. 
(TIF) 
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